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Abstract 



We derive the spin-wave dynamics of a magnetic material from the time- 
dependent spin density functional theory in the linear response regime. The 
equation of motion for the magnetization includes, besides the static spin 
stiffness, a "Berry curvature" correction and a damping term. A gradient 
expansion scheme based on the homogeneous spin-polarized electron gas is 
proposed for the latter two quantities, and the first few coefficients of the 
expansion are calculated to second order in the Coulomb interaction. 
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The study of the ground-state properties of magnetic materials within the framework 
of spin-density functional theory (SDFT) is by now a mature field By comparison, 

the study of excited-state properties is still in its infancy. There has been great interest 
in recent years in deriving a closed equation of motion for the magnetization starting from 
a first-principle description of the electrons as itinerant particles PHH], rather than from 



the time-honored Heisenberg model of interacting local moments |8HlO|. An alternative 
approach is to calculate the spectrum of spin excitations from the imaginary part of the linear 
spin-spin response function [[U]]. Our objective in this paper is to unify the two approaches 



within the framework of the time- dependent SDFT. We emphasize new aspects of the physics 
beyond the adiabatic approximation (namely, dissipation) as well as a practically workable 
computational scheme. 

Elementary spin excitations in itinerant-electron magnets fall into two groups (i) Stoner 
excitations - in which a single electron quasiparticle is spin-reversed (ii) Spin waves. 
Both types of excitations can be computed from the linear spin-spin response function 
Xij(r,r';uj) = i/3 2 J °° dte lujt {[Si(r, t), Sj(r')]), (S'j(r) is the the i component of the spin- 
density operator, with (3 = ge/(2mc) and h — 1) which determines the magnetization 
rrii(r,uj) = —(3(Si(r,u)) induced by an external magnetic field B(r', uf) at a frequency to: 

m frrt = Y.jx»M")W-,uW. (i) 

3 

Stoner excitations are distributed along branch cuts of this response function, while collective 
modes show up as isolated poles in the complex frequency plane. 



A time-dependent SDFT is ideally suited for calculating Xij ||12|| . As a first step in this 



direction one solves the static Kohn-Sham equation [13|, whose eigenfunctions and eigen- 
values determine the exact equilibrium density and magnetization. One then constructs the 
linear spin-spin response function XKS,ij( r , of the Kohn-Sham (KS) system. Because 



this is a stationary noninteracting system, the calculation can be carried out exactly |14 



Finally, the response of the physical system is calculated as the response of the KS system 
to an effective time-dependent field ~B e ff which includes, in addition to the external field 
B, a many-body "exchange-correlation" field B xc . When these ideas are cast in formulas 
one obtains the well-known mathematical relation between the matrix inverses of the exact 
response function and the KS response function, namely 

[X _1 ]ij(r, r'; u) = [x/r S _1 ]y (r, r'; u) - f xc ,ij( r , r '; ( 2 ) 
where f XCj ij(r,r';u) is the tensor that connects B a . c to the induced magnetization 

B XCti (r,u) = J2j fxc,ij{r,r';u)mj(r';uj)dr'. (3) 
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Collective spin excitations can then be obtained from the solution of the eigenvalue problem 

J [x-%<r,r>)m,(r>)dr' = (4) 

3 

where m,j(r,u) is the magnetization profile in the spin wave. The problem is to find the 
frequencies uj for which this equation has nonvanishing solutions. Because \ nas both real 



2 



and imaginary parts these eigenfrequencies will be complex in general and will determine 
both the dispersion (Re u) and the linewidth (Im uj) of spin waves. 

To appreciate the power of Eqs. (0) and (HI) we now use them to "derive" both the 
adiabatic spin dynamics [§-§1 and the Landau-Lifshitz (LL) equation [pjl in the linear re- 



sponse regime. First of all, we choose to focus on the transverse part of the response function, 
namely, the part that describes the response to a magnetic field perpendicular to the ground- 
state magnetization. It turns out that in a collinear magnet the transverse response function 
is rigorously decoupled from the longitudinal one in the absence of spin-orbit interactions. 
In general, this decoupling is justified by the difference between the time scales of the lon- 
gitudinal and transverse spin dynamics. Our key assumption is that both xlcs an d fxc can 
be Taylor-expanded, at low frequency, in powers of u>: this is justified because in a magnetic 
material the scale of the frequency dependence of x~ks is se ^ by the "Stoner-gap" , while f xc 
is controlled by multiple electron-hole pair excitations, whose spectral density is smooth. 
Keeping only the first order term in the frequency expansion of x _1 we come to 

lx~%(r, r'; u) = a^r, r') + iw^r, r') (5) 

where ay(r, r') = [x _1 ]ij(r, r'; uo = 0) is the (symmetric) spin stiffness tensor, given by the 
second derivative of the ground-state energy with respect to the magnetization, and 

(r, r ) = km , (6) 

with the derivative taken along the real frequency axis. (f2y is purely real because the 
derivative of Rex vanishes at uo = 0.) Notice that, by definition, 

aij (ry) mj (ryr' = (7) 

where E[m] is the ground-state energy regarded as a functional of m(r) The tensor 



flij can be split into antisymmetric and symmetric components as follows 

%(r,r') =%(r,r')+7^(r,r'), (8) 

where flij(r, r') = — Qji(r', r) and 7y(r, r') = 7jj(r', r). Therefore, substituting (|) into (|]), 
and switching to a real time representation with the substitution — iu — > d/dt, we obtain 
the equation of motion 

This equation reduces to the Niu-Kleinman adiabatic equation of motion 0] if the sym- 
metric tensor 7^ is neglected and the antisymmetric tensor Qij is identified with the "Berry 
curvature" . Indeed, after a lengthy but straightforward calculation we can show that 



o t <\ or I dl Pl m ) 



drrij (r') ( 
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where ip[m] is the ground-state wave function regarded as a functional of m [JIB]. Thus, 
the antisymmetric part of Eq. (§) is equivalent to Eq. (|i~0|) : the former is, however, more 
amenable to approximation and computation. 

The symmetric part of Eq. (Q), 7^, is responsible for dissipation as one can immediately 
verify by calculating the rate of entropy production at temperature T: 

dS f 5E[m] dm(r,t) 



dt J (5m(r, t) dt 

L dr dr — m — ^( r ' r 



v 



dt " JK ' J dt 

(11) 



Not surprisingly, this term is absent in a purely adiabatic theory such as that of Ref . [[| . 
We now turn to the task of approximating the right hand side of Eq. (pi). A classic 



approximation scheme is provided by the gradient expansion |L7j . In this scheme one assumes 
that the two-point function f2y(r, r') is a short ranged function of the distance |r — r'|. It is 
then permissible, if the density and magnetization are slowly varying, to expand Q as 

f2y(r,r') = O ,y[n(r),m(r)]6(r - **') 

+ fi 2) y[n(r),m(r)]V r 5(r-r') • V r /+ ... (12) 

where flo,ij [ n , m ] an d &2,ij [ n , Tn] are the coefficients of q° and q 2 respectively in the small-g 
expansion of f2^ om (q) = / (l^° m (r — r')e _ * q '^ r ~ r 'dr in a homogeneous electron gas of density 
n and magnetization m. 

We are now in a position to prove that the standard LL equation [|15[ is simply the 
zero-order approximation (i.e., the local density approximation) in the gradient expansion 
for the Berry curvature. To this end, we consider a homogeneous spin-polarized electron gas 
with the same ground state density n and magnetization mo as those of the real system at 
point r. The homogeneous magnetization is maintained by an external fictitious magnetic 
field B parallel to m , 

Bo= f de(n,m) \ ^ 



dm 



m=mo 



where e(n,m) is the energy density of the homogeneous electron gas of density n and mag- 
netization m ||18|| . The transverse spin-spin response function of this system at q = is 



[X h( q = 0,u)= — (14) 

where ujq = /3B and is the two-dimensional Levi-Civita tensor with the cartesian indices 
i, j being orthogonal to the direction of m . Thus, for the homogeneous electron gas, £lij(q = 
0) = ey/(/3m ) (see Eq. (|J)) , and then, from Eqs. @ and (|1~2|), we see that the local 
density approximation takes the form 

^•^•ra^- 1 -''' (15) 
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With this approximation Eq. @ reduces to the linearized LL equation, 



dmj(T,t) 
dt 



dE[m\ 
drrij (r) 



(16) 



Notice that the dissipative 7^ is exactly zero at this order of approximation. Thus, the 
gradient expansion for 7^ begins with a second order term 



7i 3 -(r, r') = 7 2 ,y(r)V r 5(r - r') ■ V r ' + 



(17) 



This makes physical sense because a global rotation of the spins must be rigorously un- 
damped in the absence of spin-orbit interactions. 

The remainder of this paper is devoted to the calculation of the leading gradient cor- 
rections to the LL equation. After lengthy calculations, which will be described in detail 
elsewhere, we obtain the small-g expansion of the transverse spin-spin response function of 
the homogeneous spin-polarized electron gas: 



nq 
AmrriQ 13 



(2mmo) ; 



3V 



UJ 2 — UJ 2 



| 2nn 2 (2g(0) - 1) / 

3ma \ (uj — uj ) 2 ' (uj + uj )' 

| F- + (u)(5ij -iejj) + F+_(w)((Ji i + iey) 



(uo - LU y 



(18) 



where V is the volume and ao is the Bohr radius, T-f and are the kinetic energy operators 
associated with up-spin and down-spin electrons respectively, the angular brackets denote 
the ground-state or thermal ensemble average, g(0) is the pair correlation function at zero 
separation, and F + _(uj) is a four-point response functions, defined as 



UJ 



-oj 



3V 3 



k,k' 



((5 + (-k)p(k);^_(k')p(-k'))). 



(19) 



S±(k) = S x (k) ± iS y (k) are spin-density fluctuation operators, v(k) = 4ne 2 /k 2 , and p(k) 
is the density fluctuation operator. The "Zubarev product" is defined as ((A;B)) W = 
—i / °° dte luJt ([A(t) , B]) . Taking the small uj limit of this expression we obtain the coeffi- 
cients of the gradient expansion for Q and 7 as follows: 



n 



13 (mm uj ) 2 
ReF, 



m 



(T T -Ti> , 2 7 rn 2 (2< ? (0) 



+ 



(0) 
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1 dReF. 
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3maoUJo 



[OJ 



duo 



u=0. 



(20) 
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and 

$n ImF + _(uj) , . 

72 « = ; — \n lim ^—^J- (21) 

Since the long wavelength spin wave frequency in a ferromagnet is proportional to q 2 , the 
above results indicate that the gradient corrections to Q and 7 will affect the dispersion and 
damping of ferromagnetic spin waves beginning at order q 4 . 

Equations fl2Tj|) and (BT]) contain both ground-state (thermal ensemble) averages, such 
as (Tf — Tj_) and g(0) and the dynamical response function F + _(u;) which depends on the 
spectrum of excited states. The former can be calculated to a high degree of accuracy by 
variational and diffusion Monte Carlo techniques [fL9[| ; the major challenge lies therefore in 



the calculation of F + -(u). The form of Eq. flT9D suggests that we evaluate F + _(lu) to second 
order in the Coulomb interaction: this is accomplished by substituting the noninteracting 
expression for the four point response funtion ((§+(— k)p(k); S-(\sf)p(— k')})^. Even such 
an approximate calculation turns out to be very difficult, but we have been able to establish 
analytically the limiting forms of the imaginary part of for high and low frequency. 

At low frequency ( oj « E F ^ — E F ^ where E F ^ is the largest of the two Fermi energies) we 
find 

, > m 2 Tuj(uj 2 + An 2 k 2 B T 2 ) , . 

ImF + Joj) = K —— ^ — -, 22 

where k B is the Boltzmann constant and T is the temperature. The dimensionless coefficient 
T is the sum of "direct" and "exchange" terms: T = + which we report separately 
for reasons that will become clear in the following: 

{D) _J2A_ fl(3A-l)(3A-l) 

1 - A 2 2A(1 - A) 

2 l-A 2A 1-A' K J 

where A = (1 — Q l ^/{1 + C) 1//3 > an d C = ( n t — n i)/ n ls the degree of spin polarization. 
Eq. ( P2"| ) yields the damping tensor according to Eq. (|2"T|): 

72 - = ^ 2vr^ T (3^oao) 2 (fcBT)2 - ^ 

Note that the dissipation vanishes as T 2 for T — > pOf . 

To calculate the correction to the Berry curvature we also need the real part of F + -(lo). 
To this end, we make use of the Kramers-Kronig dispersion relation 

ReF + ^) = I r d.' ImF ;- {lJ,) (25) 

7T J-00 OJ — OJ 



and use for ImF + _{uj) at finite frequency the mode-decoupling approximation of Ref. |2T 



^/m((S + (-k)p(k);S_(k')p(-k')>>, 



-g x -^- I ImxnnQt, uj - oj')Imx+-(-K oj')dio' 

IT JO 

(26) 
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where Xnn(k, uS) = V _1 ((p(k); p{— k))) w is the density-density response function, and the 
factor g x = (T^ + T^)/T^ is used to include the exchange contribution and assure the 
correct behavior in the most important low- frequency limit. 

Carrying out the calculations, we find that the last two terms on the right hand side of 
Eq. (|20P are 

^-<°> = -* w p < a (27) 



and 



dReF,(u) 



3 



= -^t^Q(C), (28) 



did 

where P(C) and Q(C) can be very accurately parametrized as 

P(C) = 1.9606 - 3.5C - 1.4C 2 InC + 2.08C 2 , (29) 

and 

Q(C) = 1.18C - 0.186C 2 - 0.842C 3 

-(0.045C-1.49C 2 )hiC (30) 

Finally, we provide practical approximate expressions for the ground-state averages ap- 
pearing in Eq. (p0|) . In the absence of detailed Monte Carlo results for intermediate polar- 
ization we suggest simple interpolation formulas between the paramagnetic (£ = 0) and the 
fully spin-polarized (£ = 1) state. For the zero-separation pair-correlation function g(0), 
a reasonable interpolation is g(0) = g (0)(l — ( 2 ) where go(0) is the paramagnetic state 
value, which is given in Ref. [Q. As for (Tf — Tj_), we propose a linear interpolation for the 
correlation contribution: 

m - T ; > = (T T - T L ) - N^±(( = 1)C, (31) 

where e c is the correlation energy per particle, r s the Wigner radius, and (Tj — T\)q = 
(kF] 5 — kFi 5 )/(20Tr 2 m) : the difference of kinetic energies of the noninteracting systems. 

The generalized Landau-Lifshitz equation (||) is the central result of this paper. It 
includes both the adiabatic spin dynamics and the conventional Landau-Lifshitz equation as 
special cases. It reduces to the adiabatic spin dynamics if the damping tensor 7 is neglected. 
It further reduces to conventional Landau-Lifshitz equation if the gradient corrections to Q 
are neglected. In addition, we have developed a density-functional scheme for the systematic 
calculation of Q and 7. The analytical expressions for the damping tensor and the Berry 
curvature make the application of the new equation of motion to spin dynamics of ferro- and 
antiferromagnetic systems quite promising. If greater accuracy is required, one can revert 
to the full-fledged linear response formalism, in which x~ks is treated exactly and only f xc 
is approximated. 

We gratefully acknowledge support from NSF grant No. DMR-0074959 and Research 
Board grant No. URB-00-029 at the University of Missouri. 
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